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High flux of electron neutrinos^) is produced at nuclear power reactors through the decays of 
nuclei activated by neutron capture. Realistic simulation studies on the neutron transport and 
capture at the reactor core were performed. The production of 51 Cr and 55 Fe give rise to mono- 
energetic is c 's at Q- values of 753 keV and 231 keV and fluxes of 8.3 x 10 -4 and 3.0 x 10 -4 ^ c /fission, 
respectively. Using data from a germanium detector at the Kuo-Sheng Power Plant, we derived 
direct limits on the u c magnetic moment and the radiative lifetime of /i„ < 1.3 x 10~ 8 fiB and 
T v /m v > 0.11 s/eV at 90% confidence level (CL), respectively. Indirect bounds on Tv/mj; were also 
inferred. The f e -flux can be enhanced by loading selected isotopes to the reactor core, and the 
potential applications and achievable statistical accuracies were examined. These include accurate 
cross-section measurements, studies of mixing angle #13 and monitoring of plutonium production. 



PACS numbers: 14.60.Lm, 13.15.+g, 28.41.-i 

I. INTRODUCTION 

Results from recent neutrino experiments provide 
strong evidence for neutrino oscillations due to finite neu- 
trino masses and mixings Q, Q- Their physical origin 
and experimental consequences are not fully understood. 
Studies on neutrino properties and interactions can shed 
light on these fundamental questions and constrain the- 
oretical models necessary for the interpretation of future 
precision data. It is therefore motivated to explore al- 
ternative neutrino sources and new neutrino detection 
channels. 

The theme of this paper is to study the production 
of electron neutrinos (v c ) from nuclear power reactors. 
Fluxes derived from the "Standard Reactor" configu- 
ration were used to obtain direct limits on the neu- 
trino properties from data taken at the Kuo-Sheng Power 
Plant. The hypothetical "Loaded Reactor" scenario was 
also studied, where selected materials were inserted to 
the core to substantially enhance the ^ c -flux. The de- 
tection channels and the achievable physics potentials in 
ideal experiments were investigated. 



II. STANDARD POWER REACTOR 
A. Evaluation of Electron Neutrino Fluxes 

Production of electron anti-neutrinos(!7 c ) due to /3- 
decays of fission products at power reactors is a well- 



studied process. There are standard parametrizations 
for the reactor i7 c spectra 3] . The typical fission rate at 
the reactor core with a thermal power of Pth in GW is 
0.3 x 10 20 Pth s _1 , while an average about 6 ;7 c /fission 
are emitted. The modeling of the v c energy spectra 
above 3 MeV is consistent with measurements at the <5% 
level0, while the low energy portion is subjected to much 
bigger uncertainties |5j- In a realistically achievable set- 
ting at a location 10 m from a core with P t h=4.5 GW, 
the z7 e -flux is 6.4 x 10 13 cm _2 s _1 . 

Nuclear reactors also produce v c via (a) electron cap- 
ture or inverse beta decay of the fission products and (b) 
neutron activation on the fuel rods and the construction 
materials at the reactor core. There were unpublished 
studies 6] on the reactor j/ -fluxes from early reactor ex- 
periments, indicating that they would not contribute to 
the background in the measurements with v e . We ex- 
tended these studies with realistic simulations, and fo- 
cused on the potentials of using them as sources to study 
neutrino physics. 

Primary fission daughters are predominantly neutron- 
rich and go through f3~ -decays to reach stability. Direct 
feeding to isotopes which decay by /^-emissions or elec- 
tron capture(EC) is extremely weak, at the ~10 _8 /fission 
level pj . The leading components for the 235 U and 239 Pu 
fissions with relative contributions rf, fission yields Yf 
and branching ratio BR for f c -emissions are shown in 
Table HJl. The average v c -y\e\d per fission Y v is therefore 
Y„ = rf • Yf • BR. In addition, stable fission products 
can undergo (11,7) capture to unstable states which de- 
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TABLE I: The leading i/ c -yields per fission (Y„ = r f • Y f • BR) 
from (a) direct feeding of daughters (Z,N) and (b) neutron 
capture on stable isotope (Z,N-1) at equilibrium conditions. 

[a] 



Series 


(Z,N) 


Y f (Z,N) Q(MeV) 


BR(%) 


Y v 




8b Rb 


1.4e" 


5 


0.53 


0.005 


4.3e" lu 


(r f = 0.62) 


87 Sr 


<le~ 


-5 


0.2 


0.3 


<1.9e _s 




104 Rh 


7e" 
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1.15 


0.45 


2.0e" 10 




128j 


1.2e" 


-8 


1.26 


6 


4.3e" 10 


23y Pu 


128j 


1.7e" 


-6 


1.26 


6 


2.6e -8 


(r f = 0.26) 


110 Ag 


1.3e~ 
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0.88 


0.3 


1.0e _s 


(b) 


Series 


(Z,N) Y f (A-l) 


cr n7 (b) Q(MeV) BR(%) Y v 




1U4 Rh 


3.2e" 2 
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1.15 


0.45 


9.0e" b 


(r f = 0.62) 


128j 


1.2e -3 


6.2 


1.26 


6 


4.3e -5 




122 Sb 


1.2e" 4 


6.2 


1.62 


2.2 


1.6e- 6 




110 Ag 


3e" 4 


89 


0.88 


0.3 


5.6e -7 


^y pu 


128]- 


5.2e -3 


6.2 


1.26 


6 


8.2e _b 


Of = 0.26) 


104 Rh 


6.8e" 2 
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1.15 


0.45 


8.0e -5 


110 Ag 


Lie" 2 


89 


0.88 


0.3 


8.9e- 6 




122 Sb 


4.3e" 4 


6.2 


1.62 


2.2 


2.6e~ 6 



cay by ^-emissions. The equilibrium yield of the major 
components^] are shown in Table^D. The leading contri- 
bution is from 103 Rh(n, 7) 104 Rh, where the yield summed 
over all four fissile isotopes is Y v = 2.1 x 10~ 4 i/ c /fission. 
However, under realistic settings in reactor operation, the 
time to reach equilibrium is of the order of 10 years, such 
that the contribution to ^-emissions by this channel is 
also small (Y v ~ 1CT 5 ). 

A complete "MCNP" neutron transport simulation^ 
was performed to study the effect of neutron capture on 
the reactor core materials, which include the fuel ele- 
ments, cooling water, control rod structures and con- 
struction materials. While the layout is generic for most 
nuclear power reactors, the exact dimensions and ma- 
terial compositions were derived from the Pth=2.9 GW 
Core#l of the Kuo-Sheng(KS) Nuclear Power Station in 
Taiwan, where a neutrino laboratorvjlOj has been built. 
The reactor core materials and their mass compositions 
are summarized in Table Ink . A homogeneous distribu- 
tion of these materials inside a stainless steel containment 
vessel of inner radius 225 cm, height 2750 cm, and thick- 
ness 22 cm was adopted. This approximation is com- 
monly used and has been demonstrated to be valid in 
reactor design studies @. Standard parametrizations of 
the "Watt" fission neutron spectra la 111 were adopted as 
input: 

n oc exp(— E/a) sinh(v / bE) (1) 

where (a,b) depend on the fission elements. The emitted 
neutron spectra for the fissile isotopes are depicted in 
Figure ^ There are on average 2.5 neutrons generated 
per fission with energy distribution peaked at ~1 MeV. 



TABLE II: The compositions for (a) the construction ma- 



terials inside the reactor core and of the containment vessel 
used in the neutron capture studies, and (b) the three isotopes 
responsible for i/ e -emissions. 

(a) 


Functions 


Materials 


Weight (kg) 


U0 2 Fuel Elements: 


Total 

235tt 
u 

238 tj 
239p u 

241 Pu 


110000 

1 V7f\ 
lO l u 

98688 
431 
84 


Non-Fuel Materials inside 

I /~itiF din m riri+ ■ 

Fuel Container 
Cooling 

Control Rod Assembly . 


Total 
Zr-Alloy 
Water 
B 4 C 

Stainless Steel 


125000 
67500 
42500 
479 
14100 


Containment Vessel: 


Stainless Steel 


910000 


T at complete insertion. 


(b) 




Materials 


Compositions (%) 
50 Cr 54 Fe 58 Ni 


Stainless Steel SUS304 
Zr-2 alloy 


0.95 4.2 
0.005 0.006 


6.3 
0.034 



The neutrons are scattered in the core and eventually 
absorbed by either the (n, fission) processes with the fuel 
elements or the (11,7) or other interactions with the core 
materials. 
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FIG. 1: The energy spectra for emitted neutrons from the 
fissile isotopes 235 U, 238 U and 239 Pu. The spectra for 241 Pu 
is approximated to be that from 239 Pu. 

An important constraint is that an equilibrium chain 
reaction must be sustained to provide stable power gener- 
ation. This is achieved by regulating the fraction of the 
control rod assembly(£) inserted into the fuel bundles. 
This constraint is parametrized by K e ff defined as the 
ratio of neutron-induced fission to starting fission rates. 
The variation of K e ff versus £ is depicted in Figure |21 
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TABLE III: The neutron capture yields Y n 


of the major 


channels at K c ff=l. The average number of neutron emitted 


per fission is EY n =2.5. 






Channel Isotope 


Weight (kg 


;) Y n 


(n,fission) on 


98688 


0.057 


fuel element 235 U 


1376 


0.62 


239 pu 


431 


0.26 


241 Pu 


84 


0.068 


EY n (fi 


ssion) = 


1.0 


^n,7J ax u 


QOCOO 

yoOoo 


u.oy 


Pnrp Rpcinn AA/fitpT" 


42519 


0.25 


10-D 

1) 




U.ZO 


50 Cr 


9.0 


u . uuuu / 


54 Fe 


26.6 


0001 8 


58 Ni 


57.6 


0.0010 


(n, 7 ) at bu Cr 


8650 


0.00016 


Stainless Steel 54 Fe 


38200 


0.00012 


Containment Vessel 58 Ni 


57300 


0.00026 


Other capture channels: 






[mainly (n,7) on other isotopes] 




0.37 


External to Containment Vessel 




0.009 


EY n (total) = 


2.5 



TABLE IV: The v c sources and their yields Y n , Y„(both in 
10 -4 ) at the reactor core. 



The equilibrium conditions require K e ff=1.0, and the dis- 
tributions of the per-fission neutron capture yield(Y n ) 
are given in Table ITTll Only 0.35% of the neutrons es- 
cape from the containment vessel, justifying that detailed 
treatment exterior to the vessel is not necessary. 
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FIG. 2: The variations of the K c jf parameter and the neutrino 
yield Y„ from 51 Cr, as functions of control rod fraction £. 



Stainless steel "SUS304" and "Zr-2 alloy" are the typ- 
ical construction materials at reactor cores, used in the 
containment vessel and control rod assembly, as well as 
in fuel element containers, respectively. These materi- 
als contain 50 Cr, 54 Fe and 58 Ni with compositions given 
in Table ITTb. Upon activation by the (11,7) reactions, 
these isotopes produce 51 Cr, 55 Fe and 59 Ni that will sub- 
sequently decay by EC and ^-emissions. Their prop- 
erties (isotopic abundance IA, (n.,7) cross-sections er n7 , 



Isotope IA(%) 


<7 n7 (b) 


Tl 


Q(keV) BR(%) 


Y n 


Y„ 


lra Rh 


4.6 t 


146 


41.8 s 


1145 


0.45 


30* 


0.14* 


50 Cr 


4.35 


15.8 


27.7 d 


753 


100 


8.3 


8.3 


54 Fe 


5.85 


2.3 


2.73 y 


231 


100 


3.0 


3.0 


58 Ni 


68.1 


4.6 


7.6e 4 y 


1073 


100 


13.0 





averaged over 18 months reactor period 



half- life ri, EC Q-value and branching ratio BR) and 
z/ -yield (Y iy =Y n -BR) are given in Table ITVl The varia- 
tion of Y„ in 51 Cr with £ is displayed in Figure [21 The 
half-life of 59 Ni is too long and thus not relevant for v e - 
emissions. The dominant reactor v c sources are therefore 
51 Cr and 55 Fe, with total yields of Y„ = 8.3 x 10~ 4 and 
3.0 x 10 -4 ^e/hssion, implying ^ c -fluxes of 7.5 x 10 16 s _1 
and 2.7 x 10 16 s _1 , respectively, at a 2.9 GW reactor. 
The total strength corresponds to a 2.7 MCi source. 

To demonstrate the validity of the simulation proce- 
dures and results, a series of cross-checks were made. As 
depicted in Figure the control rod fraction is £=8% 
at critical condition K ff=l. The relative fission yields 
of the four fissile elements are given in Table IIIII The 
neutron energy spectrum averaged over the reactor core 
volume is depicted in Figure |3] The integrated flux is 
7.6 x 10 13 cm- 2 s-\ with 26%, 52% and 22% in the 
thermal (<1 eV), epithermal (1 eV to 1 MeV) and fast 
(>1 MeV) ranges, respectively. The maximal flux at the 
center of the reactor core is about 2.5 times the average 
value. 

Comparisons were made between these results with 
industry-standard calculations and actual reactor oper- 
ation data. Agreement to within 10% was achieved. 
In particular, the important neutron capture process 
238 U(n, 7) leads to the accumulation of 239 Pu via (3- 
decays of 239 U. The yield of Y n =0.59 per fission agrees 
with the results from an independent study [Tl| . These 
consistency requirements posed constraints to possible 
systematic effects. The leading uncertainties are ex- 
pected to arise from the modeling of the reactor core 
compositions, and were estimated to be <20%. 

The process 238 U(n, 7) generates two i7 e 's from (3- 
decays of 239 U. Adding to the 6 !7 /fission from 
the fission fragments, the total ^ c -yield is therefore 
Y p =7.2 i7 c /fission, such that Y^/Yp - 1.6 x 10" 4 . In 
particular, the v -e to z7 c -e event rate ratio at the elec- 
tron recoil energy range of 300 to 750 keV is ~ 2 x 10~ 4 , 
too small to account for the factor of two excess over 
the Standard Model values in the measured i7 e -e rates 
recently reported by the MUNU experiment fl^] . 
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FIG. 3: Energy spectrum of neutron at the reactor core, 
derived from MCNP simulations. 



B. Studies of Intrinsic Neutrino Properties 

A high-purity germanium detector has collected data 
with a trigger threshold of 5 keV at a distance of 28 m 
from the core at KS Plant. Background at the range of 
1/ (kg-keV-day) was achieved [l3j. comparable with those 
from underground Dark Matter experiments. These 
unique low energy data provide an opportunity to study 
directly the possible anomalous effects from reactor v c . 
Previous reactor experiments were sensitive only to pro- 
cesses above the MeV range. While the sensitivities are 
not competitive to those from reactor ^ e 0,0| , the stud- 
ies provide direct probes on the v e properties without 
assuming CPT invariance, and cover possible anomalous 
matter effects which may differentiate v e from v e . 

The anomalous coupling of neutrinos with photons 
are consequences of finite neutrino masses and elec- 
tromagnetic form factors [li). The manifestations in- 
clude neutrino magnetic moments {Hv) and radia- 
tive decays (r„). The searches of fi u are usually 
performed in neutrino-electron scattering experiments 
U; t + c~ — > vi, + e~ . Both diagonal and transition mo- 
ments are allowed, corresponding to the cases where 
h = h and l\ ^ 1%, respectively. The experimental ob- 
servable is the kinetic energy of the recoil electrons (T) . A 
finite neutrino magnetic moment (fii), usually expressed 
in units of the Bohr magneton 



Mb 



2m e 



c = Aira 



(2) 



will contribute to a differential cross-section term given 
by||: 



f— ) 



7raLwy-T/E 
vai T 



(3) 



where a cm is the fine-structure constant, E„ is the neu- 
trino energy and the natural unit with h—c—1 is adopted. 



The quantity ixi is an effective parameter which can be 
expressed as|l5|: 



j k 



(4) 



where U is the mixing matrix and ^jk are the coupling 
constants between the mass eigenstates Vj and with 
the photon. Experimental signatures of /x; from reactor 
neutrino experiments are therefore an excess of events 
between reactor ON/OFF periods with an 1/T distribu- 
tion. 

The 18- mo nth reactor cycle suggests that the optimal 
z/ e 's are from 51 Cr, where the half-life is ri=27.7 days. 
The equilibrium flux at 28 m is 7.3 x 10 s cm _2 s _1 . 
With the actual reactor OFF period denoted by t=0 to 
t=67 days, the background-measuring OFF* period was 
taken to be from t=30 to 101 days, during which the aver- 
age residual ^ e -flux was 37% of the steady-state ON-level. 
The ON periods included data prior to reactor OFF and 
starting from t=101 days. A total of 3458/1445 hours of 
data from the ON/ OFF* periods were used in the anal- 
ysis reported in this article. 

The focus in the /z„-search was on the T=10-100 keV 
range for the enhanced signal rates and robustness in the 
control of systematic uncertainties. The v c -e scattering 
rates due to /i„ at the sensitivity level being explored 
are much larger (factor of 20 at 10 keV) than the Stan- 
dard Model rates from v e: such that the uncertainties in 
the irreducible background can be neglected^- Similar 
event selection and analysis procedures as Ref. were 
adopted. Neutrino-induced events inside the Ge target 
would manifest as "lone-events" uncorrelated with the 
cosmic-ray veto panels and the Nal(Tl) anti-Compton 
scintillators. Additional pulse shape analysis further sup- 
pressed background due to electronic noise and the de- 
layed "cascade" events. No excess of lone-events was ob- 
served in the ON— OFF* residual spectrum. A limit of 



/'// 



< 1.3 x 10" 



Mb 



at 90% confidence level(CL) was derived. The residual 
plot and the best-fit regions are displayed in Figure 

The neutrino-photon couplings probed by v-e scat- 
terings can also give rise to neutrino radiative decays: 
v \ ~~ * + 7 between mass eigenstates fj and with 
masses mj and m^, respectively. The decay rates Tjk 
and half-lives Tjk are related to jj,-^ via|lo1| 



1 

T jk 



Tjk = 



l»jk ( m ? - m k) 3 

87r m? 



(5) 



Results from oscillation experiments^, 0] indicate 
that v c is predominantly a linear combination of 
mass eigenstates v-y and vi with mixing angle dn 
given by sin 2 6>i2 ~ 0.27. The mass differences be- 
tween the mass eigenstates are Am 2 2 ~8x 10~ 5 eV 2 
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FIG. 4: Residual plots for neutrino (a) magnetic moment and 
(b) radiative decay searches with the reactor 51 Cr ^-source. 



and Ani2 3 ~2x 10~ 3 eV 2 . Both "normal" {nor.: 
ni3 ni2 > mi) and "inverted" (inv.: ni2 > mi 3> 1113) 
mass hierarchies are allowed. The v\ — > ^3 and v 2 — > z/3 
decays are allowed only in the inverted mass hierarchy 
while V2 — > v\ is possible in both hierarchies. Adopting 
these as input, the \x v limit can be translated via Eq. 
to indirect bounds of 



na 



T23 



(inv. : v x — > v 3 ) > 1 x 10 6 s/eV' 
(inv. :vi—>v z ) > 4 x 10 22 s/eV' 



T21 



(nor. + inv. : v 2 — > v\) > 6 x 10 26 s/eV 3 



at 90% CL. These limits are sensitive to the bare 
neutrino-photon couplings and are therefore valid for 
neutrino radiative decays of in vacuum. 

It is also of interest to perform a direct search of 
Ve ~~ * Vx + 7 the signature of which is a step-function 
convoluted with detector efficiencies where the end-point 
is at E v = 753 keV for u e 's from 51 Cr[T3|. As shown 
in Figure no excess of uncorrelated lone-events was 
observed in the residual spectrum from the ON— OFF* 
data. A limit of 

t v /to. v > 0.11 s/eV 



for v c at 90% CL was derived. This implies 

— > 0.08 s/eV 
mi 

— > 0.03 s/eV 
n%2 

in the mass eigenstate basis. These direct radiative de- 
cay limits apply to all the kinematically allowed decay 
channels and cover possible anomalous neutrino radia- 
tive decay mechanisms in matter, since the decay vertices 
are within the active detector volume. In particular, the 
matter-induced radiative decay rates can be enhanced 
by a huge factorf~10 23 )|l9| in the minimally-extended 
model. 

Previous accelerator experiments provided the other 
direct "laboratory" limits on the v c magnetic moments 
and radiative decay rates: a u < 1.1 x 10 -9 /ieH^] an d 
T„/m„ > 6.4 s/eV|llI, both at 90% CL. The new lim- 
its from reactor v c are complementary to these more 
stringent results, since they probe parameter space 
with lower neutrino energy and denser target density 
which may favor anomalous matter effects. Astrophys- 
ical arguments ji^ placed bounds which are orders of 
magnitude stronger |]J, but there are model dependence 
and implicit assumptions on the neutrino properties 
involved|l4). Limits were also derived from solar neu- 
trinos, through the absence of spectral distortion in the 
Super-Kamiokande spectra: u u < 1.1 x 10~ 10 fl5ll23j . 
and the observational limits of solar X- and 7-rays: 
T v /m v > 7 x 10 9 s/eV 0, both at 90% CL. However, 
the compositions of the mass eigenstates being probed 
are different from those due to v e flavor eigenstate at the 
production site studied by the reactor and accelerator- 
based experiments, such that the interpretations of the 
limits are not identical. 

The limits for the radiative decay lifetimes for mass 
eigenstates v\ and v 2 from reactor and solar neutrino 
experiments are summarized and depicted in Figure 
using the latest results from the neutrino experiments 2, 
UtI as input. The notations are defined in the figure 
caption. Several characteristic features can be identified. 
The solar neutrino experiments lead to tighter limits than 
those from reactor u e 's. The indirect bounds inferred 
from v c -e scatterings are much more stringent than the 
direct approaches, but only apply to decays in vacuum. 
Among the various approaches, only the direct limit with 
reactor v c reported in this article covers decays in both 
vacuum and matter. 



III. LOADED POWER REACTOR 
A. Enhancement of Neutrino Flux 

Using the simulation software discussed above, we in- 
vestigated the merits of inserting selected materials to the 
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TABLE V: The (11,7) and f e -yields for selected materials 
loaded to the reactor core, at K e fr=l and £=0%. 



1 1 1 1 1 1 

0.5 0.5 1 
log,,, [ m„ (eV) ] 



FIG. 5: Summary of the results on neutrino radiative life- 
times for v\ and V2 from reactor i/ e and solar neutrinos ex- 
periments, denoted by r and s, respectively. The superscripts 
(I,D) correspond to indirect bounds and direct limits, while 
the subscript "12" is attributed to decays driven by Am^, 
and so on. The upper bound (m up ) on m„ is due to limits 
from direct mass measurements, while the lower bounds m 1 ^ 
and mj° v are valid for the normal and inverted hierarchies, 
respectively. The indirect bounds 1-23 and s 2 3 are valid for 
inverted hierarchy only, while r\ 2 and s* 2 apply to v 2 — > v\ 
decays in both hierarchies. All modes are valid for decays in 
vacuum, while r D applies also for decays in matter. Bounds 
for v\ and V2 can be represented by the same bands in this 
scale. 



reactor core to enhance the i/ -flux. A convenient proce- 
dure is to load them to the unfilled rods or to replace part 
of the UO2 fuel elements or control rod assembly during 
reactor outage. Though such a scenario involves diffi- 
culties with reactor operation regulations and requires 
further radiation safety studies, it is nevertheless tech- 
nically feasible and ready — and costs much less than 
the various accelerator neutrino factories projects, which 
involve conventional neutrino beam upgrades [2^, muon 
storage rings [26) and beta beams |27j. It is therefore of 
interest to explore the physics potentials and achievable 
sensitivities. 

The candidate isotopes are those with good IA, cr n7 
and BR as well as convenient lifetimes for the activated 
nuclei. In order to sustain the fission chain reactions 
(that is, having K e fj=l), the control rod fraction £ in the 
core should be reduced and there is a maximum amount 
of the neutron-absorbing materials that can be inserted. 
This amount as a fraction of the fuel-element mass is 
denoted by A. Several selected materials and their max- 
imal A, Y n and Y„ at K e ff=l and £=0 in both natural(n) 
and pure(p) IA form are given in Table E| The optimal 
choice is 50 Cr(p). To illustrate how the allowed amount 
of control rods and source materials would relate to the 
reactor operation, the variations of K c ff and Y„ versus A 
are plotted in Figures[U<&b, at two configurations where 
the control rod fractions are (a) £=4% and (b) £=0%, re- 
spectively. Criticality condition K c ff=l requires a maxi- 
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Tl 
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Y„ 




(%) 
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(keV) 
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0U Cr(n) 


4.35 


15.8 


27.7 d 


753 


100 


14.3 


0.056 


0.056 


5 °Cr(p) 


100 










5.4 


0.31 


0.31 


M Cu(n) 


69.2 


4.5 


12.7 h 


1675 


61 


16.3 


0.20 


0.12 


63 Cu(p) 


100 










14.8 


0.25 


0.15 


lbI Eu(n) 


47.8 


2800 


9.3 h 


1920 


27 


0.073 


0.092 


0.025 


151 Eu(p) 


100 










0.035 


0.095 


0.027 



mum load of 50 Cr(p) corresponding to A=5.4% when the 
control rods are completely retrieved (£=0%). This gives 
rise to a neutrino yield of Y iy =0.31 z^/fission, and there- 
fore a Y^/Yp ratio of 0.04. As shown in Table |D] the 
total weight of non-fuel materials inside the containment 
vessel is 1.14 times that of the fuel elements. Therefore, 
such loading of 50 Cr(p) is only a small addition of mate- 
rials to the reactor core. In the case of a P t h=4.5 GW 
reactor, this maximal loading implies 8900 kg of 50 Cr(p). 
A total of 4.2 x 10 19 of i/ e 's per second are emitted from 
the core, equivalent to the activity of a 1.1 GCi source. 
The influx at 10 m is 3.3 x 10 12 cm _2 s _1 . 



a) 
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<e« 1 
0.98 ■ 
0.96 • 
0.94 
0.92 
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< 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



0.15 



0.1 



- 0.05 



A(%) 



0.3 



0.95 • 




0.2 
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FIG. 6: The variation of the K e g parameter and the neutrino 
yield Y„ from 50 Cr(p) as a function of loading fraction A, at 
control rod fraction (a) £=4% and (b) £=0%. 
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TABLE VI: Expected ^ C NCC rates per ton-year at a reactor 
51 Cr ^ c -flux of 3.3x 10 12 cm _2 s _1 (Rcore), at the standard solar 
model 7 Be Hux(Rq), and due to a 1 MCi 51 Cr source(R src ). 



Target 


IA(%) 


Threshold(kcV) 


Rcorc 


R© 


R f 


7i Ga 


39.9 


236 


2100 


3.8 


58 


100 Mo 


9.63 


168 


2300 


3.9 


64 


115 In 


95.7 


118 


11000 


19 


300 


i76 Yb 


12.7 


301 


3700 


7.4 


100 



for four half-lives of data taking. 



B. Detection and Potential Applications 

In order to detect such neutrinos, detection mech- 
anisms common to both v c and i7 c such as neutrino- 
electron scatterings are not appropriate. Instead, flavor- 
specific charged-current interactions (z^NCC) would be 
ideal. Solar u e has been observed by ^ C NCC in radio- 
chemical experiments on 37 CI and 71 Ga, with calibra- 
tion measurements using 51 Cr ^-sources performed for 
71 Ga[2^| . Detection of the low energy solar neutrinos has 
been a central topic in neutrino physics. There are many 
detection schemes and intense research program towards 
counter experiments with f c NCC[29|, using isotopes such 
as 100 Mo, 115 In, 176 Yb. The i/ c NCC rates for the various 
isotopes in their natural abundance at a 51 Cr z^-flux of 
3.3 x 10 12 cm _2 s _1 are summarized in Table fVTI Also 
listed for comparison are the rates from the standard so- 
lar model 7 Be z/ c -flux and from a 1 MCi 51 Cr source in- 
side a spherical detector of 1 m diameter. More than 10 4 
z^oNCC events or 1% statistical accuracy can be achieved 
by one ton- year of data with an indium target. Calcu- 
lations of the z/ c -flux depends on the amount of loaded 
materials and the well-modeled reactor neutron spectra, 
so that a few % uncertainties should be possible. Similar 
accuracies can be expected on the z/ e NCC cross-section 
measurements. This would provide important calibration 
data to complement the solar neutrino program. 

Such mono-energetic ^-sources and the detection 
schemes may find applications in other areas of neutrino 
physics. We outline two of such applications and derive 
their achievable statistical accuracies. Discussions on the 
systematic uncertainties and background of actual exper- 
iments are beyond the scope of this work, and will largely 
depend on the results of the ongoing research efforts to 
develop realistic ^ c NCC-detectors. 

The first potential application is on the study of the 
mixing angle #13. The mono-energetic is^s allow sim- 
ple counting experiments to be performed. The rates 
between NEAR and FAR detectors can be compared 
to look for possible deviations from 1/L 2 , L being the 
core-detector distance. The zy c -flux can be accurately 
measured by the NEAR detectors, and the oscillation 
amplitude is precisely known at fixed Am 2 , E„ and 
L. Therefore, reactor v a experiments are expected to 



TABLE VII: Sensitivities to #13 from maximally-loaded re- 
actor core with Cr(p) sources for different detector options. 
Listed are event rates per 500-ton-year(R5oo) for the FAR de- 
tector at Lo =340 m, their achievable la statistical((75oo) and 
sin 2 2#i3(<5(sin 2 2^13)) accuracies. 



Target 


IA(%) 


R500 


a 500 (%) 


<5(sin 2 20i 3 ) 


luu Mo(n) 


9.63 


1900 


2.3 


0.027 


100 Mo(p) 


100 


20000 


0.71 


0.0083 


115 In(n) 


95.7 


9100 


1.1 


0.012 


176 Yb(n) 


12.7 


3100 


1.8 


0.021 


176 Yb(p) 


100 


24000 


0.64 


0.0075 



have better systematic control than those with fission 
zZ,'s|30| where, because of the continuous energy distri- 
bution, the energy dependence in the detector response 
and the oscillation effects have to be taken into account. 
Considering both oscillation and luminosity effects, the 
sensitivities at a given neutrino energy depend on 
[sin 2 (^f^)]/%/L 2 . The optimal distance for the FAR 
detector at Am 2 =0.002 eV 2 and £^=747 keV for the 
51 Cr-source is therefore Lo =340 m. Table IVTT1 shows 
the achievable sensitivities to sin 2 2^13 with various detec- 
tor options in both natural(n) and pure(p) IA located at 
Lo from a two-core power plant each with P t h=4.5 GW. 
The source strength of Y„=0.31 f G /fission for maximally- 
loaded 51 Cr in Table El is adopted. It can be seen 
that a ~1% sensitivity can be statistically achieved with 
5 years of data taking using a 100 ton indium target — 
a level comparable with those of the other reactor- and 
accelerator-based projects. 



5 0- .. 

-0.05 - T 
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-0.15- I 

-0.2 - 

-0.25 - t 

' I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 

FIG. 7: Simulated correlations between the fractional 
changes of i/ c -yields (AY„) and those of 238 U(n,7) 239 U rates 
(AY n7 ) in the case for a 51 Cr source. Conditions under which 
the error bars are assigned are explained in the text. 

Another possibility is on the monitoring of unwar- 
ranted plutonium production during reactor operation — 
an issue of paramount importance in the control of nu- 
clear proliferation|3l) . Plutonium is primarily produced 
by /3-decays following 238 U(n, 7) 239 U whose cross-section 
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is overwhelmed at high energy(>l eV)@. In contrast, 
the (n, 7) processes in Table which give rise to v c - 
emissions are predominantly thermal. The core neu- 
tron spectra can be modified without affecting the fis- 
sion rates through optimizations of the control rod and 
cooling water fractions, making excessive plutonium pro- 
duction undetectable by monitoring the thermal power 
output alone. Measurements of the time- variations of the 
f NCC event rates are effective means to probe changes 
in the neutron spectra, and therefore to monitor directly 
the 239 Pu accumulation rates. Illustrated in Figure [7| 
are the correlations between the fractional changes of 
the ^-yields (AY,,) and those of the 238 U(n,7) 239 U rates 
(AY n7 ) in the case for a 51 Cr source having a strength 
of Y„=0.31 ^ /fission. The uncertainties in AY„ corre- 
spond to those statistically achieved with 19 days of data 
using a 10-ton indium detector located at 10 m from the 
reactor core. Such a measurement is adequate to make a 
3cr detection on a 4% reduction of the f c -flux, which cor- 
responds to a 10% enhancement of the 239 Pu production 
rate. 
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